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Abstract 

The simple measure of complexity T a p of Shiner, Davison and 
Landsberg (SDL) and the statistical one C, according to Lopez-Ruiz, 
Mancini and Calbet (LMC) , are compared in atoms as functions of the 
atomic number Z. Shell effects i.e. local minima at the closed shells 
atoms are observed, as well as certain qualitative trends of T a ^(Z) 
and C(Z). If we impose the condition that L and C behave similarly 
as functions of Z, then we can conclude that complexity increases with 
Z and for atoms the strength of disorder is a ~ and order is f3 ~ 4. 

1 Introduction 

There are various measures of complexity in the literature. A quantitative 
measure of complexity is useful to estimate the ability of a variety of physical 
or biological systems for organization. According to jlj a complex world is 
interesting because it is highly structured. Some of the proposed measures of 
complexity are difficult to compute, although they are intuitively attractive, 
e.g. the algorithmic complexity [2J |3j defined as the length of the shortest 
possible program necessary to reproduce a given object. The fact that a 
given program is indeed the shortest one, is hard to prove. In contrast, 
there is a class of definitions of complexity, which can be calculated easily 
i.e. the simple measure of complexity T a s according to Shiner, Davison, 
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Landsberg (SDL) [3J, and the statistical measure of complexity C, defined 
by Lopez-Ruiz, Mancini, Calbet (LMC) 0EIEI- 

Whereas the perfect measure of complexity is as yet unknown, the present 
work reports analysis of electron densities at atoms. We specially refer to 
the shell structure (periodicity), using two of the simplest of complexity 
measures, SDL and LMC calculated as functions of the atomic number Z. 
This is a continuation of our previous work jSj. Those measures have been 
criticized in the literature E3 IH] and a discussion is presented in Section 

m 

Our calculations are facilitated by our previous experience and results for 
the information entropy in various quantum systems (nuclei, atoms, atomic 
clusters and correlated atoms in a trap-bosons) jH], [T21 ITTA ITU Hoi fTTjl IT71 
I18[ I19[ I20[ I21j . A remarkable result is the universal property for the in- 
formation entropy S = a + b lniV where N is the number of particles of 
the quantum system and a, b are constants dependent on the system under 
consideration 13 . In fact, if one has a physical model yielding probabilities 
which describe a system, then one can use them to find S and consequently 
calculate the complexity of the system (in our case the atom) as function 
of Z. This was done in [S], where we calculated the Shannon information 
entropies in position-space (S r ) and momentum-space (Sk) and their sum 
S = S r + Sj. as functions of the atomic number Z (2 < Z < 54) in atoms. 
Roothaan-Hartree-Fock electron wave functions (RHF), for 1 < Z < 54, 
were employed |22| . For Z = 1 there is no electron-electron effect. Higher 
values of Z [Z > 54), due to the relativistic effects, are not considered. 
Analytic wavefunctions are available in |22| only for Z up to 54. Their 
importance lays in the fact that we can derive accurate wavefunctions in 
momentum space; in order to assure that the integration is accurate. In [H] 
we calculated the SDL measure with RHF densities for 1 < Z < 54. In the 
present work we calculate LMC measure in the same region of Z, for the 
sake of comparison. One could consider the case of Hartree-Fock wavefunc- 
tions for the extended region 1 < Z < 102 |23j . It turns out that both cases 
are satisfactory (work in progress). 

In [Hj and the present work, we examine if complexity of an atom is 
an increasing or decreasing function of Z. This is related with the ques- 
tion whether physical or biological systems are able to organize themselves, 
without the intervention of an external factor, which is a hot subject in the 
community of scientists interested in complexity. 

However, special attention should be paid with respect to the meaning 
of complexity or structure, which may depend on the system under con- 
sideration. In the case of atoms, the overall electrons through pair-wise 
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electron-electron interaction under the external nuclear potential, lead to 
the characteristic electron probability distribution. The selected information 
measure which shows up in this manner is employed to estimate complexity. 
Periodicity is clearly revealed here. 

In Section |21 we review briefly measures of information content and com- 
plexity. Section |21 contains our numerical results for the LMC measure and 
discussion. In Section |1] we comment on the validity of SDL and LMC 
measures. 

2 Measures of information content and complexity 
of a system 

The class of measures of complexity considered in the present work have 
two main features i.e. they are easily computable and are based on previous 
knowledge of information entropy S. 

The Shannon information entropy in position space Sr is 



where p(r) is the electron density distribution normalized to one. The cor- 
responding information entropy Sk in momentum space is defined as 



where n(k) is the momentum density distribution normalized to one. 
The total information entropy is 



and it is invariant to uniform scaling of coordinates, i.e. does not depend 
on the units used to measure r and k, while the individual S r and Sk do 
depend [TSj . 

S represents the information content of the quantum system (in bits if 
the base of the logarithm is 2 or nats if the logarithm is natural). For a 
discrete probability distribution {pi}i=i,2,...,k, one defines instead of S, the 
corresponding quantities H and -ff max 








(3) 
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and 

H max = logk. (5) 

The uniform (equiprobable) probability distribution p\ = P2 = . ■ ■ = p k = r , 
gives the maximum entropy of the system. It is noted that the value of H max 
can be lowered if there is a constraint on the probabilities {pi}. 

Another measure of the information content of a quantum system is the 
concept of information energy E defined by Onicescu [23] , who tried to define 
a finer measure of dispersion distribution than that of Shannon information 
entropy. Onicescu's measure is discussed in 

For a discrete probability distribution (pi,P2, ■ ■ ■ ,Pfe), E is defined as 

E = J2pI (6) 

1=1 

while for a continuous one p(x) is defined by 

E= I p 2 (x)dx. (7) 



One can define a measure for information content analogous to Shannon's 
S by the relation 

= i. (8) 

For three dimensional spherically symmetric density distributions p(r) and 
n(k), in position- and momentum-spaces respectively, one has 



POO 

E r = p 2 (r)Airr 2 dr, 
Jo 



(9) 



poo 

E k = n 2 (k)Airk 2 dk. (10) 
J o 

The product E r ■ E k is dimensionless and can be considered as a measure 
of dispersion or concentration of a quantum system. S and E are reciprocal. 
Thus we can redefine O as ^ 

E r Ek ' ^ ^ 

in order to be able to compare S and E. 

Landsberg j^j defined the order parameter f2 (or disorder A) as 

fi = l-A = l--^-, (12) 

^max 
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where S is the information entropy (actual) of the system and Smax the 
maximum entropy accessible to the system. It is noted that = 1 cor- 
responds to perfect order and predictability while Q = means complete 
disorder and randomness. 

In [3] a measure of complexity T a a was defined of the form 

r Qj/3 = A a nP = a q (i - a) 13 = nP(i - n) a , (13) 

which is called the "simple complexity of disorder strength a and order 
strength 0" . One has a measure of category I if /3 = and a > 0, where 
complexity is an increasing function of disorder, while category II is when 
a > 0, (5 > and category III when a = and (3 > where complexity is 
an increasing function of order. In category II complexity vanishes at zero 
order and zero disorder and has a maximum of 

(r ^ W = (a + /3)(^) atA = —p andQ = ^Tp- (U) 

Several cases for both a and (5 non-negative are shown in Fig. 2 of Ref. 0], 
where T a a is shown as a function of A. In our previous work [S] we obtained 
A = Sj Smax or $7 = 1 — A as a function of Z and plotted the dependence 
of r Qj ^ on the atomic number Z. 

We employed 5 max according to rigorous inequalities holding to atoms 
|201 [H], which hold for other systems as well (nuclei, atomic clusters and 
atoms in a trap-bosons) as verified in ^IJEZl- These inequalities are 



(15) 
(16) 
(17) 

The lower and the upper limits can be written, for density distributions 
normalized to one 



SVmin = ^ (1 + ln7r) - ^ In (t T \ 

Srmax = \ (1 + In tt) + | In Q (r 2 ) ) , (18) 
S kmi n = ^(l+lnvr)-^lnQ(r 2 )). 



»SV min 


< S r 


< 


Sr maxj 


$k min 


<s k 


< 


•Sfc maxi 


'S'min 


< s 


< 


'S'max- 



5 



S fcmax = ~(l+ln7r) + ~ ml -Tj. (!<)) 



3 ■ ' '.2v 



Smin = 3(1 +1I17T), 

S max = 3 (1 + In tt) + | In U^TI, (20) 

where (r 2 ) is the mean square radius and T is the kinetic energy. We employ 
in l|12|) S'max according to relation (f2~Ul) . 

As an alternative, we may use instead of T a p the following statistical 
measure of complexity C due to Lopez-Ruiz, Mancini and Calbet [Sj defined 
as 

C = S-D, (21) 

where S denotes the information content stored in the system (in our case 
the information entropy sum S = S r + Sk) and D is the disequilibrium of the 
system i.e. the distance from its actual state to equilibrium [7] defined 
for a discrete probability distribution {p{\ as follows jB] 



N f \\ 2 - 

D({Pi}) =^2[Pi-jf) , where Pi > and Z)w = !• ( 22 ) 
j=l ^ ' »=1 



L> is the quadratic distance of the actual probability distribution {pi} 
to equiprobability. In the continuous case, the rectangular function p[x) = 
2j;, where — L < x < L, is the natural extension of the equiprobability 
distribution of the discrete case. Thus the disequilibrium could be defined 
as 



(23) 



If we redefine D omitting the constant adding term in D* (which is very 
small for large L), the disequilibrium reads now 

D(p(x)) = J\(xfdx, (24) 

where D is positive for every distribution and minimal for the rectangular 
function which represents the equipartition. For large values of L (L — ► oo) 
relation (|24[) gives 

D = p 2 {x)dx. (25) 
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Another derivation of the formula (|25|) for continuous p(x) can be found 
in Section 3 of [2S] by using the Renyi generalized entropy 



where q is an index running over all the integer values. 

According to jSJIHlIZl S and D are the two basic ingredients for calculating 
complexity. In our 3-dimensional case, we employ instead of (|25jl the formula 



where E r , E^ are defined in Q, (|1U|) . Relation (|27|) extends the definition 
of measure of disequilibrium of the system according to LMC, to our case, 
where we are interested jointly in position- and momentum-spaces. It turns 
out that LMC definition of disequilibrium D (|25|) is identical to Onicescu's 
formula (J7J) for the information energy E. In fact, inspired by the work of 
Lopez-Ruiz, Mancini and Calbet, a new interpretation of Onicescu infor- 
mation energy may be proposed i.e. it represents the disequilibrium of the 
system or distance from equilibrium. Additionally in our case something 
new is introduced, that is the effect of a delicate balance between conjugate 
spaces, reflected in the sum S = S r + Sk and the product D = E r ■ E^. Both 
S and D are dimensionless. 

3 Numerical results and discussion 

The dependence of T a « on Z for atoms has been calculated recently in jS] . In 
the present letter we calculate C(Z) employing the same RHF wave functions 
in the same region 1 < Z < 54, for the sake of comparison. Our results are 
shown in Fig. ^and Fig. |2 We compare them with T a p(Z) shown in Fig. 3 
of [H| for (a, 13) = (1, 1), (1, 1/4), (1/4, 0), (0, 4). In all (six) cases we observe 
that the measures of complexity show local minima at closed shells atoms, 
namely for Z=10 (Ne), 18 (Ar), 36 (Kr). The physical meaning of that 
behavior is that the electron density for those atoms is the most compact 
one compared to neighboring atoms. The local maxima can be interpreted 
as being far from the most compact distribution (low ionization potential 
systems). This does not contradict common sense and satisfies our intuition. 
There are also local minima for Z=24 (Cr), 29 (Cu), 42 (Mo). Those minima 
are due to a specific change of the arrangement of electrons in shells. For 
example, going from Z=2A (Cr) with electron configuration [Ar]4s 1 3d 5 to 




(26) 



D — E r ■ Efc 



(27) 
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the next atom Z=25 (Mn), with configuration [Ar]4s 2 3d 5 , it is seen that 
one electron is added in an s-orbital (highest). The situation is similar for 
Z=29 (Cu) and Z=42 (Mo). The local minimum for Z = 46 (Pd) is due 
to the fact that Pd has a 4d 10 electron configuration with extra stability of 
electron density. It has no 5s electron, unlike the neighboring atoms. There 
are also fluctuations of the complexity measures within particular subshells. 
This behavior can be understood in terms of screening effects within the 
subshell. The question naturally arises if the values of complexity correlate 
with properties of atoms in the periodic table. An example is the correlation 
of Onicescu information content O with the ionization potential (Fig. 4 of 
[H]). A more detailed/systematic study is needed, which is beyond the scope 
of the present report. 

Our calculations in [Sj show a dependence of complexity on the indices 
of disorder a and order (3. In ||, we made a general comment that there 
are fluctuations of complexity around an average value and atoms cannot 
grow in complexity as Z increases. The second part of our comment needs 
to be modified. Various values of (a,/3) lead to different trends of T a ^(Z) 
i.e. increasing, decreasing or approximately constant. In addition, in the 
present Letter we compare C(Z) with T a ^(Z) and we find a significant 
overall similarity between the curves Tq^(Z) and C(Z) by plotting C(Z) 
and 200 x T$^(Z) in the same Fig. ©. The numerical values are different 
but a high degree of similarity is obvious by simple inspection. There is also 
the same succession of local maxima and minima at the same values of Z. 
Less striking similarities are observed for other values of (a, (5) as well, e.g. 



Concluding, the behavior of SDL complexity depends on the values of 
the parameters a and (3. The statistical measure LMC displays an increas- 
ing trend as Z increases. An effort to connect the aforementioned measures, 
implies that LMC measure corresponds to SDL when the magnitude of dis- 
order a ~ and of order (5 ~ 4. In other words, if one insists that SDL 
and LMC behave similarly as functions of Z, then we can conclude that 
complexity shows an overall increasing behavior with Z. Their correlation 
gives for atoms the strength of disorder a ~ and order (3 ~ 4. 

A final comment seems appropriate: An analytical comparison of the 
similarity of Y a ^{Z) and C(Z) is not trivial. Combining equations (|12j) and 
(|13(l we find for the SDL measure 



ri,i(Z) and C(Z). 




(28) 



S 



while for the LMC one has 



C = S ■ (E r Ek) ■ 



(29) 



S and S, 



max 



depend on Z as follows 



S(Z) = S r (Z) + S k {Z) = 6.257 + 1.069 ]nZ 



(30) 



(almost exact fitted expression) (5^1 E], while 



5max(-Z) — 5", 



r max 



(Z) + 5 fcmax (Z) = 7.335 + 1.658 InZ 



(31) 



(a rough approximation). We mention that S r , Sk, E r , E^ are known but 
different functional of p(r) and n(k) according to relations (^Q), © and ©, 
(jlOj) respectively. It is noted that our numerical calculations were carried 
out with exact values of S r , Sk, E r , E^, while our fitted expressions for 
S(Z), Smax(Z) are presented in order to help the reader to appreciate 
approximately the trend of T a ^(Z) and C(Z). 

4 Comments on the validity of SDL and LMC 
measures 

Complexity is a multi- faceted and context dependent concept, difficult to be 
quantified. It is extremely difficult at present to define complexity under all 
possible circumstances. Instead, one can choose a pragmatic approach and 
use as a starting point definitions of complexity found in the literature and 
attempt to evolve the existing framework. Thus, we have chosen the simple 
SDL measure and the statistical LMC one, which are relatively easy to 
compute. The question naturally arises if the above measures do represent 
the concept of complexity as expected from semantics, intuition or other 
general criteria. SDL and LMC measures have also been criticized in several 
aspects El ECQ , some of them are the following: 

1. All systems with the same disorder A have the same i.e. the 
function r Qj/ g(A) is universal. 

2. Landsberg's definition of disorder A = S/S m£LX does not describe prop- 
erly and does not capture the system's structure, pattern, organiza- 
tion, or symmetries. 

3. The calculation in @1 of Tn for equilibrium Ising systems is not clear. 
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In our previous and present work we use ground state electronic densities 
of atoms to evaluate T a g(Z) and C(Z). At present it is not possible to 
answer semantic or sophisticated questions as the ones raised in 1101 lllj 
and described above. One cannot claim that SDL measure is the measure 
of complexity but it is a simple measure, indicating that it can be used as a 
starting point. To our knowledge, there are no other studies of complexity 
in quantum systems as functions of the number of particles. 

The definition of r Q ( g is based on Landsberg's disorder A in terms of 
information measures *S, S mS/ ^ available as functions of TV. By definition Q 
and A are related by + A = 1. An alternative approach where A and 
O should be considered independent of each other, is an open one. Before 
considering such an extension, we calculate T at p(N) given in (|25|) . where 
S and S'max show different non-trivial dependence on N. In this sense the 
universality of T a ^(A) is modified: T a> p(N) depends on the quantum sys- 
tem under consideration, i.e. atoms |Hj, nuclei, atomic clusters e.t.c. (work 
in progress). Instead of varying the number of particles N, one could keep 
it constant, and study the effect on complexity of other parameters of the 
system. A welcome property of a definition of complexity supporting its va- 
lidity might be the following: If one complicates the system by varying some 
of its parameters, and this leads to an increase of the adopted measure of 
complexity, then one could argue that this measure describes the complexity 
of the system properly. 

There have been several criticisms regarding the LMC measure of com- 
plexity in the sense that it does not deserve the adjective statistical, it is 
not a general measure that quantifies structure and it is not an extensive 
quantity Analogous arguments with SDL could be given to support the 
use of LMC as a first working prototype of complexity to be modified by 
more sophisticated models. In various properties of the SDL measure 
are examined and a modification is attempted to remove its insensitivity 
to system differences. However in our present (quantum) case of atoms, by 
employing solely electron density distributions, it is difficult to check these 
ideas. 

In addition, the similarity of the qualitative behavior of T a r(Z) and 
C(Z) for atoms used as a case study (although they obey different defini- 
tions) is interesting. It shows that both measures share some common traits 
and correlate strongly with the periodicity of the elements in the periodic 
table. An important issue is the question if T and C are true measures 
of complexity (structure, pattern etc) or are just functionals of electron 
densities. Further research is needed to clarify their region of validity and 
applicability as measures of complexity. 
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Information entropy is an extensive quantity. The question arises if T a p 
and C are extensive or intensive quantities, what is their thermodynamic 
limit, their behavior at the boundaries (large or small order), their properties 
under scaling transformations e.t.c. For example, it is stated in ^U] that 
C is neither an intensive nor an extensive thermodynamic variable and it 
vanishes exponentially in the thermodynamic limit, for all one-dimensional, 
finite-range systems. However, in a broader context, T a ^ and C can serve as 
indicators of complexity, based on a probabilistic description. In some cases 
it can be applied in a setting wider or different than thermodynamics. Such 
a case, is the calculation of C for various shapes of probability distributions 
p{x) in 0, i.e. the rectangular, the isosceles-triangle, the Gaussian and the 
exponential probability distributions which are classified according to the 
corresponding values of C. 

Another case is the present work, where r Q;( g and C are estimated for the 
ground state electron densities, p(r) and n(k), of a quantum system (atom) 
at zero temperature. To our knowledge, such an application is done for the 
first time for a quantum system, and deserves more detailed investigation 
(e.g. for excited states). It is true that there are certain shortcomings of 
r Q ( g and C for various cases in statistical mechanics. In quantum systems 
further investigation is needed. Our title ". . .Atoms as a testbed" describes 
well our aim. 

5 Acknowledgments 

The work of K. Ch. Chatzisavvas was supported by Herakleitos Research 
Scholarships (21866) of EIIEAEK and the European Union. The work 
of Ch. C. Moustakidis was supported by Pythagoras II Research project 
(80861) of EIIEAEK and the European Union. 



11 




10 20 30 



Figure 1: The LMC measure of complexity C as function of the atomic 
number Z of atoms. 




Figure 2: Comparison of C(Z) and 200xr , 4 (Z) 
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